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Abstract 

The relationship (resemblance and/or contrast) between quantum and classical in- 
tegrability in Ruijsenaars-Schneider systems, which are one parameter deformation 
of Calogero-Moser systems, is addressed. Many remarkable properties of classical 
Calogero and Sutherland systems (based on any root system) at equilibrium are re- 
ported in a previous paper (Corrigan-Sasaki). For example, the minimum energies, 
frequencies of small oscillations and the eigenvalues of Lax pair matrices at equilib- 
rium are all "integer valued" . In this paper we report that similar features and results 
hold for the Ruijsenaars-Schneider type of integrable systems based on the classical 
root systems. 



1 Introduction 

If a many-body dynamical system is (Liouville) integrable at both classical and quantum 
levels, the classical system appears to share many 'quantum' features. For example, the 
frequencies of small oscillations near equilibrium are 'quantised' together with the eigen- 
values of the associated Lax matrices at the equilibrium. This phenomenon of close rela- 
tion/contrast between quantum and classical integrability has been explored extensively by 



Corrigan- Sasaki [ X j (this paper will be referred to as I hereafter) for Calogero-Moser (C-M) 
systems 13 El based on any root system [3J Ej- in this paper we show that similar 
phenomenon occurs in Ruijsenaars-Schneider [El El EH E] systems, which are one parameter 
deformation of C-M systems. The equations determining equilibrium can be presented in 
a form similar to Bethe ansatz equations. The equilibrium positions are described as the 
zeros of certain 'deformed' classical polynomials, which tend to the Hermite, Laguerre and 
Jacobi polynomials in the C-M limit pQ. The frequencies of small oscillations at equilibrium 
are either integers or a sort of 'deformed'' integers according to the potential is rational or 
trigonometric, respectively. In the C-M limits, these frequencies tend to those presented in 
I, which provide non-trivial support for the current results. 

This paper is organised as follows. In section 2, the Ruijsenaars-Schneider (R-S) systems 
[S] are briefly introduced and the formulas for evaluating the frequencies of small oscillations 
at equilibrium are derived. In most cases, the evaluation of the frequencies is done numer- 
ically. In section 3, the rational Ruijsenaars-Schneider systems with two types of confining 
potentials for the classical root systems A, B, C, BC and D are recapitulated. The 
formulas of the frequencies of the small oscillations are presented and compared with the 
corresponding results in Calogero systems {2] given in I. In section 4 the trigonometric R-S 
systems for the classical root systems are reviewed briefly. The formulas of the frequen- 
cies of small oscillations at equilibrium are presented and compared with the results of the 
Sutherland systems |3J given in I. 

2 Ruijsenaars-type systems 

The Ruijsenaars-Schneider systems are "discrete" versions of the C-M systems 01311], that 
is the momentum variables appear in the Hamiltonian not as polynomials but as exponential 
(hyperbolic) functions. In quantum theoretical setting this would mean that the operator (an 
"analytic difference operator", according to S. Ruijsenaars) e ±/3p = e Ttl3nB ^ dq causes a shift 
of the wavefunction by a finite unit f3h in the imaginary direction, i.e. ip(q) to ip(q =F ij3h). 
The parameter j3 has the dimensions of a momentum -1 and can be expressed as (3 = 1/mc, 
in which m and c are constants of the dimensions of the mass and the velocity, respectively. 
The R-S systems can also be considered as a one-parameter (c) deformation of the C-M 
systems (which correspond to the c — > oo limit) and are sometimes referred to, somehow 
misleadingly, as "relativistic" Calogero-Moser systems. See |T2] for comments on this point. 
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The dynamical variables of the classical Ruijsenaars-Schneider model are the coordinates 
{qj\j = l,...,r} and their canonically conjugate momenta {pj\j = 1, ...,r}, with the 
Poisson bracket relations: 

Uj,Pk} = S jk , Uj,qk} = {pj,Pk} = o. 

These will be denoted by vectors in W 

q =(q 1 ,..., q r ), p = . . . ,p r ), 

in which r is the number of particles and it is also the rank of the underlying root system 
A. In this paper we will discuss those models associated with the classical root systems, 
namely the A, B, C, D and BC. The fact that all the roots are neatly expressed in terms 
of the orthonormal basis of W makes formulation much simpler than those systems based 
on exceptional root systems. Throughout this paper we adopt the convention of /3 = 1 for 
convenience. 

Following Ruijsenaars-Schneider [Hj and van Diejen [Sj, let us start with the following 
Hamiltonian 

r 

H(p,q) = [2™s\ip^VM)V*{q)-(VM) + V*(q))] , (2.1) 

in which V* is the complex conjugate of Vj. The form of the function Vj = Vj(q) is determined 
by the root system A as: 

A: Vj(q) =w(q j )]Jv{q j -q k ), j = l,...,r + l, (2.2) 

Mi 

B, C, BC, LD: V 3 (q) = w(q J )l[v(q J -q k )v(q J +q k ), j = l,...,r. (2.3) 

The elementary potential functions v and w depend on the nature of interactions (rational, 
trigonometric, etc) and the root system A. 

2.1 Equilibrium position and frequencies of small oscillations 

It is easy to see that the system has a stationary solution 

p = 0, q = q, (2.4) 
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of the canonical equations of motion: 
in which q satisfies 



3=1 



Pi 



dHjp, q) 
dq 3 



dH(0,q) 



d qj 



0. 



By expanding the Hamiltonian around the stationary solution (J2.4j) . we obtain 
H(p, q) = K(p) + P(q) + higher order terms in p, 
in which the kinetic part K is quadratic in p 



K{p) = 

3=1 



■31 a 3, a 3 = \ V M)l 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



and the 'potential' P is given by 



P(q) = £ (2y/VM V*(q) - (VM + V*(q)) 

3=1 



3=1 V 7 



(2.9) 



This should be compared with the classical potential in C-M systems, Vc = Y^j=\(PW/ dqj) 2 , 
in which W is the prepotential jjj- It is obvious that the equilibrium is achieved at the 
point (s) in which all the functions Vj become real and positive: 



V j (q) = V j (q)*>0, j = l,2,...,r. 

Therefore the "minimal energy" P(q) is always in contrast to the C-M cases. 
Let us define the Hessian of the 'potential' P at equilibrium as Bj k : 

d 2 P{q) 



(2.10) 



dqjdq k 



j,k = l,...,r. 



(2.11) 



It is easy to verify that 



B 



It 



i dm 



~t ai dqj 



em 



dq k 



J^(?) = -i(V5(g)-VT(<z)), j,k=l,...,r. (2.12) 



Thus the small oscillations near the stationary point ()2.4|) are described by the effective 
quadratic Hamiltonian in p and q — q: 



1 r 

Heff(p, q) = K{p) + - 2J B jk (q - q)j(q - q) k . 



(2.13) 



3,k=l 
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In terms of a canonical transformation 

P'j = V^~jPv 9* = — i=(?-?)j> j = 1 ,---, r , ( 2 - 14 ) 

V 2a J 

the quadratic Hamiltonian reads 

Heff( P ',q') = ^E^) 2 + ^E 5 ^' ( 2 - 15 ) 

j=l j,k=l 

B' k = (W 2 ) lk , (2.16) 



in which a real symmetric 1 r x r matrix FT is defined by 

W jk = y/a~ 



(2.17) 



d qj 

The frequencies of small oscillations at the equilibrium are given simply by the eigenvalues 
of a matrix W: 

~ r) TV,. 

(2.18) 



which are relatively easy to evaluate. 

3 Ruijsenaars-Calogero systems 

Here we will discuss the discrete analogue of the Calogero systems [2], to be called Ruijsenaars- 
Calogero systems, which were introduced by van Diejen for classical root systems only PlllUj. 
(For the definition of C-M systems based on any root system and the associated Lax rep- 
resentation, etc, see (51E1E!-) The original Calogero systems |2j have the rational (l/q 2 ) 
potential plus the harmonic (q 2 ) confining potential, having two coupling constants g and 
to for the systems based on the simply- laced root systems, A and D, and three couplings uj 
and gi for the long roots and g$ for the short roots in B (C) root system. (For the rational 
Calogero systems, those based on B and C root systems are equivalent.) 

The discrete Calogero systems have two varieties (deformations), according to the number 
of independent coupling constants. The first has two (three for the non-simply-laced root 
systems) coupling constants g (gz and gg) and a which corresponds to u in the Calogero 
systems. This can be called a "minimal" discretisation of the Calogero systems. The second 

1 This property stems from the structure of the functions Vj, (|2.2|l . I|2.3|l and from the even nature of 
v'(x). 
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is introduced by van Diejen having three (four for the non-simply-laced root systems) 

coupling constants g (g^ and gs) and a, b both of which correspond to u. In this case the 
B and C systems are not equivalent in contrast to the Calogero systems. The integrability 
(classical and quantum) of the latter was discussed by van Diejen in some detail ITU] . 
Whereas, the former (the minimal discretisation) is new and its integrability has not been 
discussed to the best of our knowledge. As we will show in the next subsection, the very 
orderly spectrum of the small oscillations would give strong evidence for its integrability. 

3.1 Linear confining potential case 

Let us first write down the explicit forms of the elementary potential functions v and w. For 
the simply-laced root systems A and D the elementary potential functions are: 

A, D : v(x) = l — ig/x, w(x) = a + ix, (3.1) 

in which a and g are real coupling constants assumed to be positive. For the non-simply-laced 
root systems B (C), BC, they are: 

B, (C) : v (x) = 1 — igL/x, w(x) = (a + ix)(l — igs/x), (3.2) 
BC : v(x) = 1 - ig M /x, w(x) = (a + ix)(l - ig s /x)(l - ig L /2x), (3.3) 

in which g^, gu and gs are the independent positive coupling constants for the long, middle 
and short roots, respectively. As in the Calogero case, those based on B and C systems are 
equivalent. In all these cases the 'potential' P (12. 9|) grows linearly in \q\ as |g| — > oo. Except 
for the BC case, there are simple identities: 

J2i V M + VMl = const. (3.4) 

j 

Thus the Hamiltonian ()2.1|) can be replaced by a simpler one 

r 

H'(p, q) = 2j2 cosh Pj y/vMVfW, (3.5) 
i=i 

which is usually used as a starting point for the trigonometric (hyperbolic) interaction theory, 
see section HJ To be more precise, the identities are: 

r+l 

A r : i V M + V 'MY} = 2(r + l)fl + r(r + l)g, (3.6) 

i=i 
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and 

r 

B r (C r ) kD T : Yl ^) + V MT) = 2r(a+(r- l)g L + g s ) . (3.7) 
i=i 

(For D r , g L -> g and g s = 0.) 

It is interesting to note that the equations determining the equilibrium (|2.1()j) . in general, 
can be cast in a form which looks similar to the Bethe ansatz equation. For example, for 
the elementary potential ()3.1|) - (j3.2j) . (|2.1()jl read 

A: n ^*_i» = izSt, ; I r-l. (3.8) 

D r : TT ^-^-'^' + g*-^ = J = l,...,r, (3.9) 

^ <?j -qk + tg qj + qk + ig a + iqj 

IK : TT 7 ' ~ 7/ " //;/ ^ • '"- ■"' = " "li' 1 '-^. j 1 ,. (3.10) 

jhi ?j ~~ + l 0£ <?j +Qk + WL a + m qj - Ws 

They determine the zeros of deformed Hermite and Laguerre polynomials. For A r , let us 
define 

qj = VagVj, S = g/a, (3.11) 
and introduce a degree r + 1 polynomial 

r+l 

^(xH^nO*-^), (3-12) 

3=1 

which is a deformation of the Hermite polynomial H r +i(x). For lower n, Hn\x) are: 

#<P(a;) = # n (a;), n = 0, 1, 2, ^ 5) (x) = H 3 (x) - 4x5, 

Hf\x) = H 4 (x) - 32x 2 6 + 125, Hf\x) = H 5 (x) - 160x 3 5 + 200x5 + A8x5 2 , 

H^(x) = H 6 (x) -640x 4 5 + x 2 (7365 2 + 16805) - 2405 2 - 3605, ... (3.13) 

For B r let us define 

qj = VagIyj, <* = gs/gi- 1, S = g/a, (3.14) 

and introduce a degree r polynomial 

r 

L^(x) = (-l) r — yfi/rl, (3.15) 
3=1 
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which is a deformation of the Laguerre polynomial Lf (x). For lower n, we have 

r 

L^ a \x) = Li a) (x), n = 0,l, =4 a) (x) - - (-2 + 3x-3a + 2xa-a 2 ) , (3.16) 

4 M (x) = 4 a) (x) - - (- 18 + 45x- 13x 2 - 33« + 42xa -6x 2 a - 18a 2 + 9xa 2 -3a 3 

-125 + 22x5 - 22a5 + 24xa5 - 12 a 2 5 + 6 x a 2 5 - 2 a 3 5^j ,(3.17) 

Lf' Q) (x) = L^\x) - — ( - 144 + 504a; - 280x 2 + 34x 3 - 300a + 582xct - 190x 2 a 

+12 x 3 a - 210 a 2 + 210 x a 2 - 30 x 2 a 2 - 60 a 3 + 24 x a 3 - 6 a 4 
-264 5 + 760 x 5 - 241 x 2 5 - 550 a 5 + 950 x a 5 - 192 x 2 a 5 
-385 a 2 5 + 366 x a 2 5 - 36 x 2 a 2 5 - 110 a 3 5 + 44x a 3 5 - 11 a 4 5 
-144 5 2 + 300 x 5 2 - 300 a 5 2 + 420 x a 5 2 - 210 a 2 5 2 + 180 x a 2 5 2 

-60a 3 5 2 + 24xa 3 5 2 - 6a 4 5 2 ^) , ... (3.18) 

They are not the so-called u q- deformed" Hermite or Laguerre polynomials [13J. As in the 
Calogero systems, the D r is a special case gs = of the B r theory described by L^ ,_1 ^(x), 
which has a zero at x = for all r, see (1.4.20). 

It is remarkable that the spectrum of the small oscillations at equilibrium is completely 
independent of the coupling constant g, a, ql or gs- In other words, the spectrum is the 
topological invariant of the theory. It is solely determined by the root system. In fact, the 
spectrum is 

2(1 + ej ), j = l,...,r, (3.19) 
in which ej is the j-th exponent of the root system A. Explicitly, the spectrum is 

A: 2 x (l,2,...,r + l), (3.20) 

B(C)&BC: 2 x (2,4,...,2r), (3.21) 

D : 2 x (2,4, ...,2r-2,r). (3.22) 

(The lowest frequency of A series is due to the center of mass motion which is not described 
by the data of the root system.) The situation is essentially the same as in the Calogero 
systems, in which the frequencies of the small oscillations are proportional to the above 
formula ()3.19|) and are independent of the coupling constant (s) of the rational potential pQ. 
We strongly believe that this very orderly spectrum is a good evidence for the integrability 
of this type of models, as is the case for the Calogero systems. 
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3.2 Quadratic confining potential case 

Let us first write down the explicit forms of the elementary potential functions v and w. For 
the simply-laced root systems A and D, the elementary potential functions are: 

v(x) = 1 — ig/x, w(x) = (a + ix)(b + ix), a, b, g > 0. (3.23) 

The elementary potential functions are (a, b, g L , g M , g s > 0): 

B : v{x) = \ — ig^/x, w(x) — (a + ix)(b + ix)(l — igs/x), (3.24) 
C : v(x) = 1 - igs/x, w(x) = (a + 2ix)(b + 2ix)(l - ig L /2x), (3.25) 
BC : v(x) = 1 — igui 'x, w(x) = (a + ix)(b + ix)(l — igL/2x)(l — igs/x). (3.26) 

The C system is slightly different from the one given by van Diejen jlOj . since the latter is 
the quantum theory. In the limit ft, — > 0, 1)3.25)1 is the same as van Diejen's. Here again those 
based on B and C systems are equivalent in terms of the overall scaling of the potential and 
scaling of the coupling constants. This fact is reflected in their spectra ()3.29j) . (|3.3(J|) . In all 
these cases the 'potential' P ()2.9|) grows quadratically in \q\ as |g| — > oo. In the present case 
the identities ()3.4)1 are replaced by 

^ f -2(r + l)q 2 + const, A r 

except for the BC case. 

Again the equations determining the equilibrium ()2.10|) are Bethe ansatz-like equations 
which are expected to determine two-parameter deformation of the Hermite and Laguerre 
polynomials. These will be discussed elsewhere. The spectrum of the small oscillations at 
equilibrium has a very simple form. Explicitly, the spectrum is 

A: 2 j[a + b + g(r-(j -l)/2)], j = 1, 2, . . . , r + 1, (3.28) 

B: 4j[a + b + g s + gL(2r-j-l)], j = l,2,...,r, (3.29) 

C: 8j[a + b + g L + 2g s (2r-j-l)\, j = l,2,...,r, (3.30) 

BC : Aj [a + b + g L /2 + g s + g M (2r j = l,2,...,r, (3.31) 

D : Aj[a + b + g(2r-j-l)}, j = 1, 2, . . . , r - 1, 

and 2r(a + b + g(r- 1)). (3.32) 
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The Calogero models are obtained in the singular limit, a, b — > oo and by division by ab. 
In this limit, the above spectrum of small oscillations, ()3.28j) - (j3.32j) will be proportional to 
those of the Calogero models, i.e. (|3.19p . ()3.20|) - ()3.22j) . as expected. 

It is interesting to compare the above spectrum of small oscillations with the quantum 
energy eigenvalues. The quantum spectrum is given by van Diejen 0: 

A: E H = n j [n j +2(a + b)-l + 2g(r + l-j)], (3.33) 

i<j<r+i 

B: E H = A n j [n j + a + b-l + g s + 2g L (r-j)], (3.34) 

l<i<r 

C: E n = 8 J2 n i [nj + a + b-l + g L + 4g s (r-j)] , (3.35) 

l<j<r 

BC : E n = 4j2 n ih + a + b - 1 + 9l/2 + g s + 2g M {r - j)] , (3.36) 

l<i<r 

in which n = (ni, . . . , n r , (n r+ i)) is a set of 'quantum numbers' parametrising the eigenstates. 
They are non-increasing, non-negative integers (ni > ri2 > • • -n r > (n r+ i) > 0). In fact, the 
B and C formulas are special cases of the BC formula. However, the D formula needs yet to 
be derived. The r, (r + 1) independent 'lowest lying' modes corresponding to the quantum 
numbers 

n= (1,0,..., 0), (1,1,0,..., 0),..., (1,1,..., 1,0), (1,1,..., 1), (3.37) 

have energies, 

A: 2j[a + b + g(r-(j -l)/2)}, j = l,...,r + l, (3.38) 

B: Aj[a + b + g s + g L (2r-l-j)], j = l,...,r, (3.39) 

C: 8j[a + b + g L + 2g s (2r-l-j)], j = l,...,r, (3.40) 

BC : Aj [a + b + g L /2 + g s + g M {2r - 1 - j)] , j = 1, . . . , r, (3.41) 

which are exactly the same as (|3.28|) - (|3.31|) . 

4 Ruijsenaars-Sutherland systems 

The discrete analogue of the Sutherland systems [3], to be called Ruijsenaars-Sutherland 
systems, was introduced originally by Ruijsenaars and Schneider [H] for the A type root 
system. The quantum eigenfunctions of the A type Ruijsenaars-Sutherland systems are called 
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Macdonald polynomials jT^j, which are a one-parameter deformation of the Jack polynomials 
P2J. Here we will discuss Ruijsenaars-Sutherland systems for all the classical root systems, 
A, B, C, D and BC [TU]. The structure of the functions {Vj(q)}, ()2.2|) and (j2.3|) are the same 
as in the Ruijsenaars-Calogero systems, but the elementary potential functions v and w are 
trigonometric instead of rational. Because of the identity Y7j=i {Vjio) + Vjil)*} = const., 
f!3.4|) . the Hamiltonian (|2.1j) can be replaced by a simpler one 

r 

H'(p, q) = 2 Y, cosh Pj y/Vj^Vfiq), (4.1) 
j=l 

which is obviously positive definite. This is also valid for the BC case in contrast to the 
rational potential cases discussed in the preceding section El 
The elementary potential functions v and w are: 

A, D : v(x) = cosh 7 — i sinh 7 cot x, w(x) = 1, (4.2) 
B : v(x) = cosh 7^ — 2 sinh 7^ cot x, w(x) = cosh7 S — 2sinh7sCotx, (4.3) 
C : f(x) = cosh 7s — i sinh 75 cot x, w(x) = cosh 7^ — z sinh 7^ cot 2x, (4.4) 
C" : f(x) = cosh 7 S — i sinh 75 cot x, w(x) = (cosh 7^ — i sinh 7^ cot 2x) 2 , (4.5) 
i?C : f(x) = cosh 7 M — i sinh7 A ,/ cotx, 

w(x) = (cosh7 S — i sinh 75 cot x) (cosh 7^ — i sinh 7^ cot 2x), (4.6) 

in which 7^, jm and 75 are the positive coupling constants for the long, middle and short 
roots, respectively. Both C and C and BC cases are special cases of the most general 
integrable interactions including the long roots [o? L = 4) introduced by van Diejen |llj . 
Note that in our paper only the classical dynamics h — > is discussed. That is, van Diejen's 
constant 7 = i(3h/2 (eq.(2.7) in jTTj) is treated as vanishing. 

For the systems based on A type root system, the above identity ()3.4|) 

r+l 

A r ■ y~] v j(Q) = sinh[(r + 1)7]/ sinh 7 (4.7) 

is known in a different context For the other root systems the identity (|3.4|) reads 

{2 sinh [r7 L ] cosh [ (r - 1 ) 7 L + 75] / sinh j L , B r , 
2sinh[r 75 ]cosh[(r-l)7 5 + 7 L ]/sinh7 5 , C r , (4.8) 
sinh[(2r- 1)7]/ sinh 7 + I, D r . 
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It is easy to verify that the B r formula reduces to the D r one for 7^ — > 7 and 75 = 0. 
Had we started from the simplified Hamiltonian (j4.1|) instead of the original one ([2.1)1 . the 
above constants ()4.7|) - ()4.8|) would give the minimal energies. In the simply-laced A and D 
cases, the r.h.s of ()4.7|) and ()4.8|) can be interpreted as " 'q- deformed integer" version of the 
dimensionality of A r and D r vector representations, [r + l] q and 1 + [2r — l] q . 
The equilibrium of the A type theories is achieved at " equally-spaced" 

q = vr(0, l,...,r- l,r)/(r + 1) + £(1,1, ...,1), £ G R : arbitrary, (4.9) 

configuration, as in the original Sutherland systems. In all the other cases, the equilibrium 
and the frequencies of the small oscillations are determined by solving the Bethe ansatz-like 
equations ()2.1()|1 numerically. Certain one-parameter deformation of the Jacobi polynomials 
is expected to describe the equilibrium for B, C, D and BC systems, which will be 
discussed elsewhere. The spectrum of the small oscillations at equilibrium has a very simple 
form. Explicitly, the spectrum is 

A : 4 sinh[(r + 1 — j)*y] sinh^'7]/ SU1I17, j — 1, 2, . . . , r + 1, (4-10) 

B : 4sinh[(2r - 1 - j)7 L + 75] sinh[j7 L ]/sinh7 L , j = 1, 2, . . . , r - 1, 

2sinh[(r - 1)7^ + 75] sinh[r7 L ]/ sinh7 L , (4.11) 

C : 4sinh[(2r - 1 - j)7s + 7l] sinh[j7 5 ]/sinh7 S , j = 1, 2, . . . , r, (4.12) 

C : 4 sinh[(2r - 1 - j) ls + 2 7i ] sinh[j 7s ]/ sinh 7s , j = 1, 2, . . . , r, (4.13) 
BC : 4 sinh[(2r - 1 - j)<y M + 2 7l + 75] smh\jj M }/ sinh 7M , j = 1, 2, . . . , r, (4.14) 
D : 4 sinh[(2r — 1 — 7)7] sinh[j7]/ sinh7, j = 1, 2, . . . , r — 2, 

2sinh[(r — 1)7] sinh[r7]/ sinli7, twofold degenerate. (4.15) 

The spectrum of the A system ([4.10J1 2 reflects the symmetry of the Dynkin diagram j «-> 
r + 1 — j. The twofold degeneracy of the D spectrum (14.15(1 also reflects the symmetry of 
the D Dynkin diagram. 

The original Sutherland models are obtained in the singular limit in which all the coupling 
constant(s) become infinitesimally small: < 7, Jl,Jm,Js 1. In this limit, the spectrum 
of small oscillations at equilibrium will be linear in the coupling constant (s). In these limits, 
the known spectrum of small oscillations obtained in I is reproduced. As for A r , the spectrum 
2 This formula was known to S. Ruijsenaars |17j . 
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()4.10j) becomes 

A: 4(r + l-j)j 7 , j = l,...,r + l, (4.16) 

which is eq.(5.16) of Corrigan- Sasaki pQ, to be referred to as (1.5.16). For B r , the spectrum 
(I4.11j) in the limit is 

B r : 4[(2r-l-j)7L + 75]j, j = l,...,r-l, 2[(r - l) 7i + 7 <?]r, [2], (4.17) 

which is (1.5.74). For the spectrum ()4.13|) in the limit is 

C": 4[(2r-l-j) 7s + 2 7i ]j, 3 = ±,---,r, (4.18) 

which is (1.5.81). For D r , the spectrum ()4.15|) in the limit is 

D: 4(2r-l-j)j 7 , j = 1, ■ ■ ■ , r - 2, 2r(r - 1) 7 , [2], (4.19) 

which is (1.5.60). The limiting spectra provide non-trivial supporting evidences for the 
formulas ljOTf|) - <|07i|l . 

It is a challenge to understand the connection between the quantum spectrum of the 
Ruijsenaars-Sutherland systems, eq. (5.17) of ^H] and the above spectrum of small oscilla- 
tions ()4.1()j) - (j4.15j) . It would be interesting to evaluate the eigenvalues of the Lax matrices 
at equilibrium as shown for the C-M systems £Q. However, the knowledge of the Lax pairs 
for the Ruijsenaars-Schneider systems is still quite limited [THj . 
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